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We study lattice wave functions obtained from the SU(2)i Wess-Zumino-Witten conformal field 
theory. Following Moore and Read’s construction, the Kalmeyer-Laughlin fractional quantum Hall 
state is defined as a correlation function of primary fields. By an additional insertion of Kac-Moody 
currents, we associate a wave function to each state of the conformal field theory. These wave 
functions span the complete Hilbert space of the lattice system. On the cylinder, we study global 
properties of the lattice states analytically and correlation functions numerically using a Metropolis 
Monte Carlo method. By comparing short-range bulk correlations, numerical evidence is provided 
that the states with one current operator represent edge states in the thermodynamic limit. We 
show that the edge states with one Kac-Moody current of lowest order have a good overlap with low- 
energy excited states of a local Hamiltonian, for which the Kalmeyer-Laughlin state approximates 
the ground state. For some states, exact parent Hamiltonians are derived on the cylinder. These 
Hamiltonians are SU(2) invariant and nonlocal with up to four-body interactions. 

PACS numbers: 71.27.+a, 73.43.-f, 11.25.Hf 


I. INTRODUCTION 

The discovery of the fractional quantum Hall (FQH) 
effecfPl has revealed the physical existence of a new, 
strongly correlated state of matter. One surprising prop¬ 
erty of FQH phases is the fact that they cannot be clas¬ 
sified in terms of symmetric^. This is a fundamental 
difference to most other known states of matter and led 
to the concept of classifying states in terms of topological 
ordeil^. Topological phases are not characterized by a 
local order parameter but their properties depend on the 
topology of the system. For example, the ground state 
degeneracy of a FQH system was shown to depend on 
the genus of the surface on which the state is definecP. 

A central characteristic of topological order is the ap¬ 
pearance of gapless edge states. Such edge modes are 
known to be present in a FQH systenJ^ even though the 
bulk is gapped. It was shown that these edge states can 
be used to characterize the topological order of a FQH 
stat^^. Furthermore, they are of particular physical im¬ 
portance because they determine the transport proper¬ 
ties of the system. The strongly correlated state formed 
by the edge excitations of a FQH system is that of a chi¬ 
ral Luttinger liquicP and its theoretical properties are re¬ 
lated to Kac-Moody algebra^^. For continuous systems, 
edge state wave functions were constructed a s corr elation 
functions of a conformal field theory (CFT)P^^ and in 
terms of Jack polynomial^^^. 

The use of I + I dimensional CFT in describing FQH 
states was pioneered by Moore and ReacP^. In this 
approach, a trial wave function for the FQH state is 
given as the correlation function of conformal primary 
fields. Laughlin’s wave function for the continuunP^ 
as well as its SU(2) symmetric, b osonic lattice version, 
the Kalmeyer- Laugh lin (KL) stat^^^^, were shown to 
be of this fornP^ElHini, Jt was already conjectured by 


Moore and ReacP^ that the FQH edge modes should 
be described by the same CFT that defines the bulk 
wave function. The idea of this bulk-edge correspon¬ 
dence was later deve loped into an extension of Moore and 
Read’s methocP^^. Within this approach, trial FQH 
edge states are obtained from CFT descendant states. 

This paper is based on a similar ansatz for edge states 
in lattice systems. Our starting point is the SU(2)i 
Wess-Zu mino-W itten (WZW) theory, which was used 
previousl^ff^^ to obtain the KL state as a correlation 
function of primary fields. By insertion of Kac-Moody 
currents, we define a tower of states, which corresponds 
to the CFT descendant states. Since we work in a spin 
formulation, the Hilbert space of the lattice system is 
that of N spin-1 degrees of freedom, where occupied sites 
are represented as spin-up and empty sites as spin-down. 

We show that the mapping from CFT states to spin 
states is surjective, i.e. any state of the spin system can 
be written as a linear combination of states constructed 
from the CFT. As a consequence, not all states obtained 
in this way are edge states. For some of the wave func¬ 
tions, we carry out numerical calculations to test if they 
describe edge states: On the cylinder, we compare spin 
correlation functions in the states with one current op¬ 
erator to the K L state using a Metropolis Monte Carlo 
algor it hrrP^^^. We show that the nearest-neighbor bulk 
correlations approach each other exponentially as the 
thermodynamic limit is taken. This indicates that the 
states with one current operator indeed describe edge 
states. 

In the pa st years, parent Hamiltonians of the KL 

stat cJ^^ I ^^ -^^^ and its non-Abelian generalization^^^HUl 
were constructed. It was also shown that t he KL stat^^ 
and non-Abelian FQH lattice state^^^^ have a good 
overlap with the ground states of certain local Hamil¬ 
tonians. For one of these Hamiltonians, an implemen- 
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tation scheme in ultracold atoms in optical lattices was 
proposecP^. In this work, we study the local Hamilto¬ 
nians of Ref. [ 30 ] on the cylinder and use the KL state 
and the edge states with one current operator of order 
one as ansatz states. We find that there is a choice of 
parameters for which both the ground state and some 
low-energy excited states have a good overlap with our 
ansatz states. 

Furthermore, we construct Hamiltonians for which the 
KL state and some of the states obtained from it by in¬ 
sertion of current operators are exact ground states. The 
proposed parent Hamiltonians are SU(2) invariant and 
nonlocal with up to four-body interactions. 

This paper is structured as follows: In Sec.jlTj we intro¬ 
duce the CFT model, define the map to a spin-^ system, 
and show that it is surjective. In Sec. |ml properties of 
our states on the cylinder are studied and numerical ev¬ 
idence is given that the states with one current operator 
describe edge states. We consider a local model Hamil¬ 
tonian in Sec. jlVj and derive exact parent Hamiltonians 
for some states on the cylinder in Appendix The con¬ 
clusion is given in Sec. |V[ 

II. CFT MODEL AND SPIN-| STATES 

In this section, we review some properties of the SU(2)i 
Wess-Zumino-Witten (WZW) theory and define the cor¬ 
respondence between states of the CFT and states of a 
spin-^ system on the lattice. It is shown that this map 
of CFT states to lattice spin states is surjective. 

A. SU(2)i Wess-Zumino-Witten theory 

We consider the chiral sector of the SU(2)i WZW the¬ 
ory. In addition to conformal invariance, this theory 
exhibits an SU(2) symmetry generated by the currents 
J^{z) with a G {x^y^z}. The modes are defined in 
terms of the Laurent expansion 

oo 

r{z)= Y. 

n= —oo 

and satisfy the Kac-Moody algebrsP^ 

[j^, = iSabcJm+n (2) 

Here, Sahc is the Levi-Civita symbol and 5 ah the Kro- 
necker delta. For indices c G {x,^, z}^ we adopt the con¬ 
vention that indices occurring twice are summed over, 
unless explicitly stated otherwise. 

The SU(2)i WZW theory has two primary fields, 
the identity with conformal weight h = 0, and a two- 
component field (j)s{z) {s = ±1) with conformal weight 


h = ^. The field 0 s(2:) provides a spin-| irreducible 
representation through its operator product expansion 
(OPE) with the SU(2) current 

~ --3— Y ( 3 ) 

s'=±l 

where ^ are the SU(2) spin operators. 

The field content of the SU(2)i WZW theory can be 
represented in terms of the chiral part ip{z) of a free boson 
field as 

J^{z) = --^d^<j){z), 

J±(z) = r{z) ± ijy{z) = : . (4) 

Here, q = 0 if the operators act on a state with an even 
number of modes of the h = | primary field and g = 1 
otherwise, and the colons denote normal ordering. The 
value of s G { — 1 , 1 } equals two times the spin-z eigen¬ 
value of (j)s{z). 

A general state in the identity sector of the CFT 
Hilbert space is a linear combination of states 

( 5 ) 

where | 0 ) is the CFT vacuum and are positive integers. 
The sum k = defines the level of the state. By 

means of the Kac-Moody algebra of Eq. ([^, a basis can 
be chosen for which the mode numbers are ordered: rii > 

ni-i > • • • > ni > 0. 

B. Spin states on the lattice 

To each CET state {J-m • • • T^ni)(0)|0)5 we associate 
a state 

ICr.'.'.ni)= E (6) 

Si,...,SN 

in the Hilbert space of a system of N spin-| degrees of 
freedom. Its spin wave function is defined as the CET 
correlator 

= {<l>sAzi) ■ • ■ ^sAzN)iJX ■ ■ • (7) 

where (...) denotes the expectation value of radially or¬ 
dered operators in the CET vacuum. 

In the sum of Eq. Si = ±1 and l^i,..., Sat) is 
the tensor product of eigenstates |s0 of the spin op¬ 
erator at position i {t^\si) = yls^)). The com¬ 
plex coordinates Zi are parameters of the wave function 
, sjsf) and define a lattice of positions in 
the complex plane. Since we want to keep them fixed, 
we do not explicitly indicate the parametric dependence 
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of (<^ 1 : • • • 5 Oil fhe positions Zi for simplicity of 

notation. 

The wave function corresponding to the CFT vacuum 
is given bjisa 

• • • 5 ^n) = ( 0 S 1 4 ^SN i^n)) ( 3 ) 

N 

= S.X.Yl{zi- 

i<j 

where 6 s is 1 if Ylf=i ^ 0 otherwise. Xs denotes 

the Marshall sign factor, 

N 

9=1 

It follows from the alternating sign in the definition of 
0 s(z) [cf. Eq. Q]. Due to the presence of the Marshall 
sign factor, the state 00 is a singlet of the total spirP^ 

T 00 o) = 0 . ( 10 ) 

We assume that N is even, as required by the charge 
neutrality condition = 0 in the wave function 

of 00- If w as sho wn previously that 0 o is equivalent to 
the KL stat d^'^ l ^^ -^. We note that the wave function 0 o of 
Eq. ^ contains a square root, which has two branches. 
It is assumed that the same branch is taken for all fac¬ 
tors in Eq. Which of the two branches is chosen is 
not important, since this choice only influences the wave 
function by a total factor. 

The states t)e related to 0 o by an appli¬ 

cation of spin operatord^ In Eq. Q, one can write the 
modes as an integral according to Eq. jQ, deform 
the integral contour, and apply the OPE (IW between 
current operators and primary fields. Defining 

N 

ul ^ (11) 

it then follows that 

(12) 

We note that the operators defined in Eq. 0 sat- 

isfy 

=^eabc«m+„, (13) 

which is a Kac-Moody algebra with vanishing central ex¬ 
tension. 

In addition to the states 0 ^ 5 consider the 
wave functions obtained by the insertion of two addi¬ 
tional primary fields, one at zq = 0 and one at Zoo = oc: 

= {4>s^ ioo)4>s^ (^l) • • • 4>sm izN)4>so (0)) 

N N 

OC (5sXs(-1)(1-®~)/2 ^s^os„/2 J-J 

n=l n<m 

( 14 ) 


where 6 s is 1 if sq + Soo + ~ ^ ^ other¬ 

wise. The state 00°’^°° contains a singlet and a triplet 
obtained by the tensor product decomposition of the two 
additional spin-^ fields. As we show in Sec. IIIB, the 
singlet component on the cylinder can be derived from 
the wave function of N primary fields on a torus in the 
limit where the torus becomes a cylinder. 

We note that the states 0 o, 0 n 5 y.‘.n 0 of 

Eqs. 0111 ). and ( p! 4 | ) are not normalized. Whenever the 
norm of a state is needed, it will be explicitly included in 
the corresponding expression. 


C. Completeness of spin states 


One may ask if the linear transformation that maps 
GET states {J-m • • • '^-ni)( 0 )| 0 ) states 0 ^ 5 ;;;^^ is 

surjective, i.e. whether any state in the Hilbert space 
l-Ljsf of N spin -1 particles can be written as a linear com¬ 
bination of the states We now show that this is 

indeed the case. Introducing the N x N matrix 



/ 1 

1 \ 




z = 

to 

. (zn)-^ 



(^iv)-(^-^V 


the definition of Eq. 0 becomes 


f < \ 


{n\ 


= z 

n 

0-(iv-i)/ 


Vn/ 


( 15 ) 


( 16 ) 


The determinant of Z is the well known Vandermonde 
determinant: 


N 


det (Z) = n(^ 


-1 




( 17 ) 


i<j 


Z is therefore nonsingular if all positions Zi are distinct, 
which we assume to be the case. As a consequence, the 
relation of Eq. ( 16 ) can be inverted, i.e. all spin operators 


tj can be written as linear combinations of with 
n G { 0 ,..., V - 1 }. 

The Hilbert space 'Hat is spanned by the states ob¬ 
tained from any nonzero state by successive application 
of spin operators Given that can be expressed 
in terms of it follows in particular that the states 
ICr.'.'.ni) = '^-ni ■ ■ ■ U-n, l^o) with m S {0, . . . , - 1} 

span 'Hat. Since 0 o is a singlet and 
total spin, a state 


( 18 ) 

can be written in terms of states for which all mode num¬ 
bers are greater than zero by commuting Uq to the right 













4 


until it annihilates [cf. the commutator of Eq. (HU)]- 
Therefore, 4 ^at is spanned by the states 


rii > 0 . 

This argument shows that not all states constructed 
by the insertion of current operators can be edge states 
compared to t/^q. For the states with one current operator 
numerical evidence will be given in Sec. IIIC that 
these represent edge states on the cylinder. 

The fact that the states span Hat raises the 

question about the minimal level k = Yl]=i needed 
to obtain the complete Hilbert space. We note that an 
upper bound is given hy k = N{N — 1), since any product 
of spin operators can be reduced to a product of at 
most N spin operators. Each of these spin operators can 
then be expanded in terms of the operators with 
nG {0 ,...,A/'—1}. We carried out numerical calculations 
for the states which indicate that the states up 

to level k = {N/2)‘^ are enough to obtain the complete 
Hilbert space from ipQ. 


III. PROPERTIES OF STATES ON THE 
CYLINDER 


In this section, we study properties of the states defined 


in Sec. HB on the cylinder. 

We consider a square lattice with lattice sites in 
the open direction and Ny lattice sites in the periodical 
direction of the cylinder. After mapping the cylinder to 
the complex plane, the coordinates assume the form 


Zn = e 




(19) 


Here, jx G {1, • • • ^Nx} is the x-component of the index 
and jy G {!,...,is the ^-component, so that j = 
{jx — 1)A^ jy ranges from 1 to A = For the 

remainder of this paper, we adopt a two-index notation, 
where this is convenient, i.e. we may write Zj^j^ instead 
of Zj, denoting the x- and ^-components by subscripts. 

Note that one has the freedom to rescale the coordi¬ 
nates since this changes the wave functions only by a total 
factor. The constant factor that we included in Eq. (19) 


is chosen such that the center of the cylinder is at the 
unit circle. 

We assume that the number of sites A is even. It is 
also possible to study the case of A being odd which will 
show the existence of two topological sectors. However, 
we can already identify the two anyonic sectors for A 


even. As will be shown below in Sec. HIB the state 
and the singlet component of the state with two ad¬ 
ditional spins (one at 2 )o = 0 and one at ^oo = c>c) can be 
obtained from the wave function of A primary fields on 
the torus in the limit where the torus becomes a cylin¬ 
der. This argumentation shows that the two states jjo 
and represent the two anyonic sectors in the case of 
an even number of spins. It would be possible to con¬ 
sider an odd number of sites on the cylinder by putting 
an additional spin either ai zq = 0 or at Zoo = oc so that 


the charge neutrality condition is satisfied. On the torus, 
however, such a construction is not possible and the ar¬ 
gumentation that we used to identify the two sectors for 
even A does not directly apply. 


A. Global transformation properties of GET states 


In this subsection, we study global transformation 
properties of the states jjo, • This 

serves two purposes: First, it allows us to conclude that 
states with a different momentum are orthogonal, i.e. 
they have different global properties. Later, we will study 
their local behavior numerically and compare spin corre¬ 
lation functions in the bulk. The symmetries derived in 
this subsection will be exploited in our numerical calcu¬ 
lations to obtain efficient Monte Carlo estimates. 

We consider the translation operator in the periodical 
direction Ty and the inversion operator X. Their precise 
definition and the derivation of their action on the states 
-00 5 are given in Appendix [ a} Geo¬ 

metrically, the translation operator rotates the system in 
the periodical direction and the inversion operator corre¬ 
sponds to a reflection of the cylinder along its two central 
cross sections. We call it an inversion because it acts on 
the coordinates defined in Eq. (19) as ^ z^^. 

Eigenstates of Ty and X are given in Table As we 
show in Appendix A 3, applying the inversion X to 
corresponds to inserting the current operators at Zo© = cxo 
instead of zq = 0. We use the notation for 

these states: 


= ( 20 ) 

Since the momentum in the periodical direction Py is re¬ 
lated to Ty through the relation Ty = we conclude 
from Table |l] that an additional insertion of a current op¬ 
erator into the correlation function of primary fields 
adds a momentum of —27Tn/Ny to the state. In particu¬ 
lar, the states jjo and have a different momentum 

if k = Yl]=i is different from 0 modulo Ny. 


B. Relation to the KL states on the torus 

In this subsection, we place the system on the torus 
and take a limit in which the torus becomes a cylinder. 
We show that the wave function of A primary fields on 
the torus gives rise to 'ipo and the singlet component 
of on the cylinder. 

We define the torus for cji > 0 and uj 2 = iL with L > 0 
by identifying a complex number 2 : with z-\-nuJi-\-muj 2 for 
m^n G Z. The two circumferences of the torus are there¬ 
fore given by uji and |cc; 2 |- Let us denote the positions on 
the torus by i.e. we assume that Vi lie in the rectan¬ 
gle spanned by uji and uj 2 . Keeping the positions Vi fixed 
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TABLE L Eigenstates of the translation operator Ty and the 
inversion T. The sum of mode numbers '^j is denoted 

by k. Eor the states , the current operators are 

inserted at Zoo = oo [cf. Eq. (|20|]. 


Eigenstate 

Ty 

T 

'00 



Yni ...n\ 


— 

M.-ai , 10.1 ...ai 

<Pni...ni V-m -ni 

— 

(±l)(-l)^>'f 

4’^ - 4’^ 

4’'^ ^4’^ + 4’’^^ 

(_l)JVxf+JV, 

(_l)JVxf+JVx 

(_l)JV«f+JVx 

+iVx + l 


and '0 q’^ 


and taking the circumference L = \uj 2 \ oo transforms 
the torus into a cylinder, as illustrated in Fig. 



FIG. 1. (Color online) Limit in which the torus becomes a 
cylinder: The circumference |ct; 2 | is taken to infinity while the 
positions lie in the finite region of size Nx x Ny (red patch). 


On the torus, there are two states with k G 

{0, These are given 


V’fe (si,...,Siv) = .■.(f>SNi'>^N))k 

/ N \ ^ _ 

X!CiSi,2T JJ 6> 2 


oc 5sXsO 


a = l 


Kj 


SiSjl2 


( 21 ) 


Here, = Vi/uJi are the rescaled positions, t = 002 !^! is 
the modular parameter of the torus, and 0 the Riemann 
theta function defined as 


0 


a 

b 


iC,r) = 

n^Z 


gi7rr(n+a)^+27ri(n+a)(C+6) 


( 22 ) 


The limit UJ 2 which transforms the torus into a 

cylinder, implies r ^ ioo. In this case, only the terms 
with the smallest value of (n + a)^ contribute to the sum 
of Eq. (22). These terms have n = 0 for a = 0 and 


n G {—1,0} for a = |. Therefore, 


(C^^r) ^ 1 , 


(C,2r) + e-^^^), and 


1 iCi - 

2 


for r ^ ioo. With follows that 

N 


n 


3 * 2 {Cm~\~Cn)SrnSn _ g* 2 '^n=l 


(23) 


(24) 


In the limit uj 2 ioo^ we therefore obtain 

V>r-(si,...,SAr) 

SmSnl‘^ 


N 

oc 4Xs II (( 


,27riVml^l _ 


m<n 


and 


t/>r-(si,...,sjv) 

2 

/ N 


oc Ssxs I n 

\n=l 

N 

-no 

m<n 


^TriVnSnll^l 


N 

n< 

n=l 


^-TviVnSnluJi 


27riVml^l _ luJi 


Sm^n 1“^ 


(25) 


(26) 


The exponentials ^ cylinder of circum¬ 
ference oji. We therefore identify Zn = and 

Ny = uji. Comparing the expressions for in the 

limit Co ’2 ^ oc to Eqs. (§ and ( p!4| ), we conclude that 

^torus ^torus ^ = ^sgl^ 

2 

It would be very interesting to investigate the relation 
between excited states in one and two dimensions on the 
circle and the plane, respectively, to those on the torus. 
The result for the relation between and is a 

2 

first step in that direction. 


C. Spin correlation functions and edge states 

We calculated two-point spin correlation functions in 
the states 7 /^ 0 , and in the singlet state 

i’f = 4’^ - 4’^ (27) 

using a Metropolis Monte Carlo algorithm. This allowed 
us to compare properties of the states numerically for 
large system sizes by sampling the relevant probabil¬ 
ity distributions. We furthermore exploited the transla¬ 
tion and inversion symmetries of Table [l| to average over 
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jx 


jx 



FIG. 2. (Color online) Two-point spin correlation function jx, ^y) in the bulk (upper panels) and at the edge (lower 

panels) for Nx = 13 and Ny = 20. The left panels show the two-dimensional dependency in a color plot. Whenever the 
value of the correlation function does not differ from zero by more than three times the estimated error, we excluded the data 
point from the plot (blank fields). In the right panels, the absolute value \S^^^{ix,ix, ^y)\ of the correlation function along the 
^-direction is plotted. Points for which the sign of the correlation function is positive (negative) are shown in blue (red). For 
the data shown in gray, the mean value does not differ from zero by more that three times the estimated error. In the bulk, 
the correlations decay exponentially, while a nonzero, negative correlation remains at the edge for Ay > 5. 


equivalent correlation functions, thus obtaining a faster 
converging Monte Carlo estimate. 

In this subsection, we shall use the notation 


are encoded in the spin correlation functions, we com¬ 
pared the nearest-neighbor two-point correlators in the 
bulk for different system sizes. The relative differences 










(V’lV’) 


(28) 


for the two-point correlation function in a state ^ip. Since 
all wave functions that we consider have a translational 
symmetry in the periodical direction, their value only 
depends on the difference Ay of the positions in the y- 
direction. 

Before comparing the wave functions with each other, 
we discuss the spin ordering pattern in - 00 , which is en¬ 
coded in the correlation function 5'^^ (i^,, A^). (Since 

ijjQ is a singlet, xx-^ yy- and 2 : 2 )-correlations are the same 
and only correlation functions with a = b are nonzero.) 
Our numerical results are shown in Fig. In the bulk 
of the system, we observe a ring-like structure with an 
alternating magnetization. At the edge, the correlations 
are anti-ferromagnetic at short distances. At larger dis¬ 
tances along the ^-direction, however, the sign becomes 
stationary and a negative correlation remains. In the 
two-dimensional picture, the ordering is still character¬ 
ized by an alternating magnetization with the sign of the 
correlation function changing along the x-direction. 

We now discuss the question if the states pJi can be 
considered as edge states. If so, then the local properties 
of tpi and ipo in the bulk should be the same. Since these 


jx^ ^y) 


^^0 {^x^jxj ^y) 


(29) 


are shown in Fig. for a = 6 = c = z (left panels) and 
a = z^b = c = X (right panels). Correlation functions for 
other choices of a, 6 , and c either vanish or can be reduced 
to these due to the SU( 2 ) invariance of pjQ. In the upper 
panels, the correlations along the ^-direction are shown 
and in the panels of the lower two rows along the x- 
direction. We find that the relative differences approach 
zero exponentially as a function of Nx. Even though the 
differences tend to be larger for smaller Ny, they are still 
exponentially suppressed as Nx is increased. This is an 
indication that the wave functions pJi indeed describe 
edge states compared to pJo as the thermodynamic limit 
in the open direction is taken. 

Our results for the comparison between and pjQ are 
shown in Fig.|^ Since both pjQ and are singlets, it is 
enough to compare the 2 ) 2 )-correlations. Furthermore, the 
correlations in the positive and the negative x-direction 
are the same in the middle of the cylinder since p jy and 


III A 


are symmetric under the inversion, cf. Sec. 

In contrast to we find that the thermodynamic limit 
in the x-direction is not enough for the differences to 
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1 5 9 13 1 5 9 13 

N, N, 


- Ny=A Ny = 12 

Ny = 8 - = 16 


FIG. 3. (Color online) Comparison of nearest-neighbor bulk 
correlations in 'ipi and 'ipo. The vertical axes show the relative 
differences \S^^z — S^^\/\S^^^\ with b = c = z (left panels) 
and h — c — X (right panels). Along the horizontal axis, the 
number of spins in the open direction (iVx) is varied. The 
different curves correspond to configurations with different 
Ny. The insets show for which sites the correlation functions 
were computed. The sites in the central column of the insets 
correspond to the middle of the cylinder in the a:-direction. 
The relative difference decreases exponentially in Nx. 


vanish. Rather, we observe that the differences become 
stationary if Ny is held fixed and increased. As shown 
in the right panels of Fig. the differences do, however, 
tend to zero exponentially as a function of Ny if N^ is 
chosen large enough. 


D. States at a higher level 


In the previous subsection, the states at level one in 
current operators were considered. We also compared 
spin correlations in to those in -00 for higher values 
of n. For very large mode numbers n, only the terms at 
the edge contribute to the sum in To see this, let 



FIG. 4. (Color online) Comparison of nearest-neighbor bulk 
correlations in and '0o- The vertical axes show the relative 
differences |/|^Jo I- panels, Nx is varied 

Wq 

along the horizontal axes and the different curves correspond 
to different choices of Ny. On the right panels, Ny varies along 
the horizontal axes and Nx is fixed. For Nx large enough, the 
differences tend to zero exponentially as a function of Ny. 


us consider 


^—n—mNy 


^ 1 

/ ^ n-\-mNy 3 

j=l ^3 


(30) 

iVx 

Ny 




■ WZ - 4 - 0 , 

Sx.jV 

Jo; — 1 

3y — ^ 

(31) 


For large values of m, the terms with > 1 are ex¬ 
ponentially suppressed with respect to those that have 
jx = 1. We denote the corresponding states with one 
current operator by Xn- 


y 

Ixn) = IC+mJvJ « ^ (32) 

3y — ^ 

Fig. ID shows the difference in nearest-neighbor corre¬ 
lations along the ^-direction relative to -00 for N^ = 13 
and Ny = 8. The three curves correspond to the states 
Xi- As the position in the open direction is in¬ 
creased, the differences vanish exponentially for all three 
states. We note that the differences are large at the left 
edge {ix = 1) and small at the right edge {ix = 13). This 
agrees with the expectation that the operators are 
localized at the left edge. In contrast to the state 
the states perturb only at one edge and their be¬ 
havior is therefore expected to approach that of at the 
other edge. The results of Fig. provide an indication 
that describe edge states also for n > 1. 















































































































E. Inner products of states from current operators 





_ xlj = xpl - 'tp = ip^ - t = Xi 


FIG. 5. (Color online) Relative difference 5^^ l/I^S ^I in 
nearest-neighbor correlators for -0 e [cf- Eq. (32) 

for the definition of Xi]- The position ix in the (r-direction is 
varied along the horizontal axis. The plots in the left panels 
have a — b — z and those in the right panels a — b — x. In 
the upper panels, the correlations along the y-direction are 
shown [jx — ix , = 1] and the lower panels correspond to 

correlations along the (r-direction [jx = Zx + 1, = 0]. 


We note, however, that the linear span of for n G 
{1,..., A/" — 1} contains not only edge states. The states 
can even be linearly combined so that is perturbed 
at an arbitrary position j : 


N 




n=2 


(33) 


where Z is the matrix defined in Eq. ( [l^ . This obser¬ 
vation can be understood from the fact that two states 
and %p^ are not necessarily orthogonal if m — n = 0 
modulo Ny. In this case, and have the same mo¬ 


mentum in the ^-direction, as discussed in Sec. |III A] The 
linear combination 




jx -1 




27r(j'a, 1) 


) E (34) 


jy = l Jx,Jy 


for example, receives no contribution from spin operators 
at the left edge {jx = 1). Even though both %p^ and 
perturbed from ipo mostly at jx = 1, this is 


not the case for the difference of Eq. (34). 


In this subsection, we discuss the relation of inner 
products between the states "ipnp.'.'^ni 
spin system and GET inner products between states 
{J-m • • • '^-ni)(0)|0)- states in the continuum 

that are constructed from descendant states of a GET, 
the authors of Ref. m come to the remarkable conclusion 
that, in the thermodynamic limit and under the assump¬ 
tion of exponentially decaying correlations in the bulk, 
the inner products between edge states are the same as 
the inner products between GET states. We now consider 
inner products between states constructed from current 
operators to test if a similar correspondence holds for the 
lattice states 'ipnp.'.'^x GET states they are con¬ 

structed from. The spin system inner products that we 
consider are given by 

d/c+/c' 

(V^olV^o) 

(^olt/'o) 

(35) 


where N is the number of spins. 


R= min \zj\ = e 


(36) 


is the minimal absolute value of the positions, k = 
and k' = The factor accounts 

for the scaling of the operators with respect to a 
rescaling of the positions. The minimal value is chosen 
because the operators 


N .a 
a _ ^j 


(37) 


have the highest contribution at the edge with \zj\ = R. 

We compare the inner products of the lattice sys¬ 
tem ( [ 3 ^ to the GET inner products 

(0|j“p..j“vX---TWl0)- (38) 


If a correspondence similar to that of Ref. |TT| also holds 
for lattice states, then the expressions of Eq. (35) should 
approach those of Eq. (38) in the thermodynamic limit. 

Note that the inner products of the spin system are 
hard to evaluate for large system sizes, whereas the GET 
inner product can be easily computed using the Kac- 
Moody algebra On the level of the spin system, the 
insertion of current operators corresponds to an applica¬ 
tion of spin operators to ipo [cf. Eq. (12)]. Therefore, 
the inner products can be determined numerically using 
a Monte Garlo method if the number of current operators 
is small. 

We calculated the inner products for the states 
and ipi I , which are all nonzero states at levels one and 
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two. For these states, inner products between different 
states vanish because they have either a different spin or 
a different momentum. It is thus sufficient to compare 
the norm squared of a state to the norm squared of the 
corresponding CFT state, as summarized in the following 
table: 

Spin state CFT state Norm squared of CFT state 
J%\0) (JoJ^o) = 1 (no sum over a) 


For a given system size, we observe a smaller differ¬ 
ence for the states at level k = 1 than for those at level 
k = 2. The computed inner products approach the CFT 
expectation if Ny is increased. The dependence on the 
number of spins in the x-direction is, however, very weak 
for Nx > 3. In particular, the CFT result is not ap¬ 
proached if Nx is increased and Ny kept fixed. For large 
enough Ny, our data suggest that the spin system inner 
products approach the CFT result with a power law in 

Ny. 


N. 





10 20 Ny Ny 


V, 


Mill I ■ ■ iM I . I I lU 

10-1 10 ° 10-1 10 ° 


10-1 10° 



FIG. 6. (Color online) Inner product of spin system states: 
Relative difference to the CFT expectation [cf. Eq. ([3^]. 
The colors correspond to the states \ (red), 'ip 2 (blue) and 
'ipi (green). The upper panels show the relative difference 
in a color plot as a function of Nx and Ny. We observe a 
very weak dependence on Nx for Nx > 3. The lower panel 
shows the dependence on Ny for Nx = 11. For large enough 
Ny, the data are consistent with a power-law behavior with 
an exponent of approximately —1.1. Monte Carlo error bars 
are not plotted because they are barely visible on the chosen 
scale. The maximal relative error of the shown data points is 
0.31 %. 


In Fig. our numerical results are shown for the rel¬ 
ative difference 

- (V’cftI^cft) 

(V^olV^o) _ 

('0CFt|'0CFt) 

as a function of the system size. Here, 'ip G {'ipi, 7/^2 ? il 
is one of the spin states, k = 1 for k = 2 for 

pj G 5 V^CFT is the CFT state corresponding 

to pj. 


IV. LOCAL MODEL HAMILTONIAN 


In the previous section, we provided numerical evi¬ 
dence that the states with one current operator insertion 
represent edge states with respect to In this section, 
we study a set of local Hamiltonians on the cylinder. For 
a suitable choice of parameters, the ground state of the 
corresponding Hamiltonian has a good overlap with ipQ 
and some of its low-energy excited states are well ap¬ 
proximated by pJi, the states with one current operator 
of order one. 

We study the local Hamiltonian^^ 

H=J2Y. + ^2 E + -^3 E 

ihj) {{id)) {id,k)o 

(40) 

In these sums, the sites lie on a square lattice, {i,j) 
denotes all nearest neighbors, ((i, j)) all next-to-nearest 
neighbors, and {i,j, k)^ all triangles of nearest neighbors 
for which i,j, and k are oriented counter-clockwise. It 
was shown in a previous stud}® that the ground state 
of H on the plane (open boundary conditions in both di¬ 
rections) and on the torus has a good overlap with the 
KL state for a range of parameters J 2 , -^25 Js- Here, 
we study on a cylinder of size Nx x Ny, where Nx de¬ 
notes the number of sites in the open direction and Ny 
the number of sites in the periodical direction. In the 
following, we parametrize H in terms of two angles 0i 
and O 2 '. 


J 2 = cos(6>i)cos {O 2 ), 

J 2 = sin (^ 1 ) cos {O 2 ), 

J3 = sin(6>2). (41) 


For Nx = 5 and Ny = 4, we studied the overlap be¬ 
tween pjQ and the ground state pjQ of as a function of 0i 
and O 2 using an exact numerical diagonalization method. 
We also computed the overlap of the states with one cur¬ 
rent operator insertion at level one %pi and the first ex¬ 
cited states of H that have spin one and the same 
momentum in the ^-direction as pjf. Here, m G { —1,0,1} 
denotes the eigenvalue of tp^. 

We denote the overlap between two states (pi and p 2 
as 




\{(t>l\4'2)\ 

WMM’ 


( 42 ) 
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p = 0 p = 1 p = 2 p = 3 



012 012 012 012 


spin 5 


FIG. 7. (Color online) Overlaps of states constructed from 
CFT and eigenstates of the local Hamiltonian H of Eq. (40) 
for Nx = b and = 4. The angles 0i and O 2 parametrize the 
coupling constants of H according to Eq ( |41[ ). In the upper 
panel, the overlap = |('0 g|^o)V(II'^g'II ll'0o||) be¬ 

tween ^0 and the ground state of H is plotted. The lower 
panel shows the overlap between 'ipl and the first excited state 
'ipE of H with the same spin and y-momentum as [spin one, 
= 0, momentum 3/(87r)]. The point marked with an open 
circle has 61 = 0.0275 x 27t and 62 = 0.06 x 27r and the highest 

combined overlap of ^Q H ~ 1.2858. 


where ||0|| = ^/{^\^ is the norm of a state. In Fig. 
the overlaps and O shown as a 

function of the parameters of the Hamiltonian. Due to 
SU(2) invariance, it is sufficient to consider the overlap 
between the states and , which both have = 0: 

i(V’ti^f)l = l(V’^^lV’r)i = i(V’ii^?)l, (43) 

where = ipf ^ The best value for the combined 

overlap ^ was obtained for the 

angles 0i = 0.0275 x 27r and O 2 = 0.06 x 27r: 

0.9829 0.8289 0.9991 0.9907 

Since the size of the Hilbert space grows exponentially 
with the system size N, the overlaps are expected to scale 
exponentially in N. By taking the Nth root, one obtains 
a measure for the overlap per site, which takes into ac¬ 
count this exponential scaling. Notice that the overlap 
per site is higher than 99 % for both the ground and the 
excited state. 

The low-energy spectrum of H for the parameters with 
the best overlaps is plotted in Fig. We find 8 en¬ 
ergies below the energy of The spectra plotted in 


FIG. 8. (Golor online) Low energy spectrum of the Hamilto¬ 
nian H of Eq. ( [40| for Nx = 5, Ny — — 0.0275 x 27r, and 

O 2 = 0.06 X 27r. The ground state energy of is Eo ~ —24.11. 
The four panels correspond to the four sectors of y-momentum 
p/{271 Ny) with pG {0,1,2,3}. Each shown level has a degen¬ 
eracy of 2s + 1 corresponding to the values of . The labels 
show the ansatz state constructed from GET and the value 
for its overlap with the corresponding eigenstate of H. The 
8 energies shown in red are those that are smaller than the 
lowest energy in the p = 1 and p = 3 sectors. [At the level 
marked with an asterisk (*), there are actually two energies 
with a splitting of approximately 1.211 x 10“^. This is not 
visible on the scale of the plot.] 


Fig. I^are separated into sectors of different p-momentum 
p/{27rNy) = p/{S7r) with p G {0,1,2,3}. Note that the 
states 2 p^ are the first excited states with p = 3. The 
spectra for the momenta for p = 1 and p = 3 are the 
same because H is invariant under the inversion opera¬ 
tor I introduced in Sec. MM 

= (44) 

The relation 




-1 


(45) 


between X and translation operator in the p-directio n Ty 
follows directly from their definition (cf. Appendix A 2 
and |A3| ). Therefore, if \^p) is an eigenstate of H with mo¬ 
mentum p/{2'KNy)^ then X\^p) is also an eigenstate with 
momentum {Ny — p)/{27rNy). This means that for 
with p = 3, there is a corresponding eigenstate X\^p^) 
with p = 1, which satisfies 




(46) 


Here, = X\^pl) is the state obtained by inserting 

the current operator at Zoo = 00 instead of zq = 0 [cf. 

Eq. @]. 
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Our results show that the state and the states with 
one current operator of order one are good approxima¬ 
tions of low-energy eigenstates of H for = 5, Ny = 4, 
0i = 0.0275 X 27 r, and O 2 = 0.06 x 2tt. This raises the 
question if further eigenstates of H are effectively de¬ 
scribed by states constructed as CFT correlators. We 
also computed the overlaps of eigenstates of H with some 
additional states constructed from current operators for 
higher orders in current operators. At level two in cur¬ 
rent operators, the overlaps with the first excited states 
that have the same spin and momentum as our ansatz 
states are given by 0.5486 for ^2 (0.9704 per site) and 
0.3301 for -0^ 5 (0.9461 per site). 

In this Section, we considered a local model and com¬ 
puted overlaps between ansatz states constructed from 
CFT and eigenstates of the model Hamiltonian. In Ap¬ 
pendix we also derive exact, SU(2) invariant parent 
Hamiltonians for some of the states constructed from cur¬ 
rent operators. These Hamiltonians are nonlocal with up 
to four-body interactions. 


V. CONCLUSION 

This work studies trial wave functions for lattice FQH 
states constructed as chiral correlators of the SU(2)i 
WZW CFT. To each CFT state, characterized by a se¬ 
quence of current operators, we associated a correspond¬ 
ing state lattice system. For continuous 

systems, analogous states constructed from CFT were 
proposed as FQH edge states previousl^ff^mi. The fact 
that we work on the lattice allowed us to apply Monte 
Carlo techniques to test a central expectation for edge 
states: That the local, bulk properties of different edge 
states should be the same. 

For a system on the cylinder, we compared spin corre¬ 
lation functions in the states with one current operator 
(V^f) to the state with no current operators (V^o)- Our nu¬ 
merical results show that the nearest-neighbor bulk cor¬ 
relations approach each other exponentially as the num¬ 
ber of spins in the open direction {N^) is increased. On 
the other hand, the states and are different glob¬ 
ally since their spin and momentum are different. This 
supports the assumption that they describe edge states. 

We compared inner products of lattice states at levels 
one and two in current operators to CFT inner products 
of the corresponding descendant states. For large enough 
Ny (periodical direction), the computed inner products 
approach the CFT expectation with a power law in Ny. 
This suggests that there is a correspondence between in¬ 
ner products of states and CFT inner products 

in the thermodynamic limit. Such a correspondence was 
found for continuous wave functions in Ref. nn 

Furthermore, we compared nearest-neighbor bulk cor¬ 
relations in to those in where is the singlet 
component of the state obtained by insertion of two ex¬ 
tra primary fields. In contrast to we find that the 
correlations do not approach each other if the thermody¬ 


namic limit is taken only in the open direction. However, 
if Nx is chosen large enough, the difference in correlation 
functions vanishes exponentially as a function of Ny. 

We showed by an exact diagonalization that has a 
good overlap with the ground state of a local Hamiltonian 
and with the first excited states that have the same 
spin and momentum as 'ip^. This could be an indication 
that further low-energy excitations of that local Hamilto¬ 
nian are edge states described by the SU(2)i WZW CFT. 
It would be interesting to investigate this relation in more 
detail for larger system sizes and different topologies. 

We showed that the complete Hilbert space is covered 
by the linear span of the states therefore, 

only a subset of these states are edge states. For the 
states with one current operator, we argued that not all 
linear combinations of states describe edge modes 
because states with the same ^-momentum can be non- 
orthogonal. It is possible to restrict the space of states to 
an orthogonal subset given by with m G { I , . . . , Ny }. 

Taking the limit of large mode numbers could be an¬ 
other possibility of removing bulk states for the linear 
span of More precisely, one can replace rii by 

Ui + rriiNy and then take rrii ^ 00 . In this limit, the 
sum in the operators only extends over the 

edge sites because all other positions are exponentially 
suppressed. The fact that this class of states (and also 
linear combinations of such states) is obtained from ipQ 
by application of edge spin operators only, suggests that 
their complete span represents edge states. For one of 
these states, we did numerical tests that indeed indi¬ 
cated that Xi is an edge state. 
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Appendix A: Translation and inversion of states on 
the cylinder 

1. Transformation under a permutation of the spins 

Both the translation operator Ty and the inversion op¬ 
erator X act on a product state as a permutation of the 
spins. Such a permutation operator Or is defined for the 
permutation r of N elements as 

Or\si, . . . , Sn) = |<Sr(l), • . • , Sr(N))- (Al) 

The action of Or on pjQ and can be rewritten 

in terms of a permutation of the positions which will 
facilitate our calculations for Ty and X. Our derivation 
of this transformation rule follows Ref. m 
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We consider the wave function 

N 

^Zl,-;ZN (^Si,...,Sn) = SsXs 11(2:* - (A2) 

i<j 

which is equivalent to 2 po because it only differs by a 
spin-independent constant [cf. Eq. ^]. We have also 
explicitly written out the parametric dependence on the 
positions Zi. Similarly, the wave function 

N 

n=l 

N 

X l[{Zn- + (A3) 

n<m 


is equivalent to . Let us first calculate the transfor¬ 
mation of under a simultaneous permutation of 

both the spins and the coordinates. Since every permu¬ 
tation can be decomposed into a series of transpositions, 
we consider the case that r is a transposition: 


r{i) = < 


h 

n, 

m, 


if i ^ {m, n}, 
if i = m, 
if i = n. 


(A4) 


for m,n G { 1 ,..., N} and m < n. It follows that 


We note that Eq. (A7) is not valid for the original wave 


function ipo but only for which differs from ipQ 

by a factor depending on Zi. However, this factor does 
not depend on the spins. Therefore, if is an 

eigenstate of O^-, then this is also the case for 

Compared to ' 00 , there are some an additional factors 
present the wave function - 0 ^o,soo, 2 i v, 2 Ar^ Since these are 
invariant under a permutation of both the spin s and the 
coordinates, a formula analogous to Eq. (A7) holds for 


'^Oj'^oo i^\ j — 
0 


2. Translation in the periodical direction 


The translation operator Ty is defined through the per¬ 
mutation Ty\ 


'Ty (i^C, iy) 


{px-i T 1 ), if iy ^ Tly^ 
{pxi f)? if 


(A8) 


where ix is the x-component and iy the ^-component of 
an index i. 

The signature of this permutation is given by 

signify) = (-l)^Aiv.-i) = (A9) 


where we used that N = NyNx is even. In terms of the 
positions, the transformation corresponds to a multipli¬ 
cation by a phase, = e‘^^^^^yzj. Therefore, 


N 

n (-1 

transformation of Vg 

m )/2(_l)( + l)/2 

n—1 

X ^_2^('®j('®rn,+Sn)+2)/2 

j=m+l 

= -1. (A5) 


= sign(7^)t^7’'*"’’ , Siv) 

i<j 

= , sn). (AlO) 

Here, we have used that 

N 

(All) 

i<j 


Therefore, if r is a general permutation corresponding to 
A/V subsequent transpositions. 


= sign(r)t^o'’-’^'^ (si,..., sjv), (A6) 


where sign(r) = (—is the signatur e of the permuta- 
tion. Substituting Si by in Eq. (A 6 ), we arrive at 

the final transformation rule: 


= sign(r)q"‘'>’'''’''"‘"’ (si,..., sjv)- (A7) 

The transformation under a permutation of the spins can 
therefore be calculated by considering the correspond¬ 
ing transformation of the coordinates and taking into ac¬ 
count the signature of the permutation. 


which follows from — 0 * The eigenvalue of ipo 

with respect to Ty is therefore (—With 

(A 12 ) 

it follows that the eigenvalue of '^nl'.'.'.ni 
^-27nk/Ny ^ where k = 

Eor we obtain 

JS0,Soo,Zi,...,Zn r 1 i 

AT , N , 

n=l n<m 

= (A13) 
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In the last equation, we have used that + <^0 + 

Sqq — 0. 


3. Inversion 


We require that the inversion X acts on the positions 
defined in Eq. (19) as 




(A14) 


This leads to the definition 


^ ij^x i^y) 


T 1 '^xi^y ^y)") ^y ^ ^y") 

{Nx ^1-ix, Ny), if iy = Ny. 


(A15) 


have 




• • • ,Sjv) 



N 

II 

0 

••’""(si,.. 




m<n 

II 

0 

••’""(si,.. 

) g ^ + O-Og^Zm Zn)~\~7ri) 

= C” 


SN)e~^ ^m,n(^rnSn-\-l){\og(ZmZn)+7Ti) 

1 

X 02 

Ei(21og(2^ 

)+7ri) 

= C” 


,. ,sjv)(-l)^^”’''''''^'’. (A19) 


In the last step, we have used that si-\ - \- sn = 0 and 


N 

E = TTi (AT + Ny;). (A20) 

n=l 

Therefore, the eigenvalue of 2 po with respect to X is 

sign(i)(-l)^^-+^- = (-1)^^^. (A21) 


We note that in our choice of Zi , the center of the cylinder 
is at the unit circle. If this is not the case, then the 


definition of Eq. (A15) leads to an additional factor when 


X is applied to Zi. 

In order to determine the sign of the permutation, we 
arrange the state |si,...,SAr) in a matrix: 


51.1 

52.1 


«2,iV, 


\ 


ki.i’ • • • ^•SiV^.iV^) - 

••• SN:^,NyJ 

The transformed state is then given by 


(A16) 


I|si,l, . . . , SN^,Ny) 

( 


SNy„Ny-l SN^,Ny-2 

SjV^-l,iVy-l SN^-l,Ny-2 ■■ 

\ Si^Ny-1 Si^N^-2 


■ SAT^,! SAfx.iVy 
• <SJV,-1,1 Sv 


Sl,l 


ShNy / 

(A17) 


To bring the transformed matrix back to the original 
form, we first reverse all Ny columns and then reverse all 
Nx rows excluding the last element of each row. A single 
sequence of L elements can be reversed in ^L{L—1) steps. 
Therefore, the sign of the permutation is given by 

sign(X) = (_ ^Ny—^)Ny—‘^) ^ (A18) 


We next determine the contribution from t he coo rdi- 
nate part of the wave function ipQ. Using Eq. (A14), we 


The states eigenstates of X, but trans¬ 

form as 




= (-l) 

Here, we have used that 

N 


""‘'<‘,•••< 1 ^ 0 ). (A22) 


7—1 — 1 , gj ^—1 _ ^ 


N 




ii) 

(A23) 


i=l 


if the center of the cylinder is at the unit circle. In terms 
of the states defined in Eq. (20), we therefore 

have 


I\rn\::Z) = (A24) 

Eor the transformed states ), the current op¬ 

erators are therefore inserted at Zoo = oc instead of at 
zq = 0. Eigenstates of X with eigenvalues (±1)(—1)^^^ 
are then given by 


(A25) 

Einally, we determine the transformation of 
^so,So^,zi,...,ZN respect to X. As for -00, there 

is a contribution from the sign of the permutation 
and from the transformation of the coordinates. The 
calculation is similar to that for '00, only that now 
So + Soo + Si = 0- We find 

— ^.*0,^1 

° (A26) 
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Since and only differ by a spin- 

independent factor, we also have 

(A27) 

Note that X exchanges the spins sq and Soo in 


Appendix B: Exact parent Hamiltonians 


As shown in Sec. IV, the edge states ipi have a good 
overlap with low-lying excited states of a local model, 
for which ipo approximates the ground state. In this sec¬ 
tion, we analytically construct SU(2)-invariant, nonlocal 
parent Hamiltonians for some linear combinations of the 
states Hamiltonians for which they are exact 

eigenstates with the lowest energy. 


1. Construction of parent Hamiltonians 


Before describing our results, we note that the condi¬ 
tion =0 for all a implies that the state -0 is part of 

the subspace on which and commute. To show 

this, we first note that 

[V^^,T’’]=iSabeK, (B4) 

whic h is a di rect consequence of the definitions of 
Eqs. ( |B1| ) and ( |B2| ). It then follows that 

[T'’T^P“] 1^) = {-iSbacT^V^ - iebacKT”) IV-) 

= -i£bac [Dn^T^] 1^) = £bac£chdDi\lp) 

= 0, (B5) 

where we assumed that = 0 for all a. The states 

satisfying = 0 can therefore be decomposed into 

sectors of different total spin. 

We note that the condition [T^T^,P^] \^Ij) = 0 is 
equivalent to 

[T^Tb^a + (I - N)T^] I'lp) = 0, (B6) 


The starting point of our construction is the operator 
N 

+ i£abct\^)- (Bl) 

i^j 

In Appendix we explicitly compute the action of 
on states constructed from ipo by insertion of current op¬ 
erators, and show that does not mix the states 

with different levels k = Yl]=i if ^ This prop¬ 

erty is key to our construction of parent Hamiltonians: 
It allows us to treat the levels separately starting with 
the lower levels, which have fewer states. The action of 
on states at level k is described by a matrix. For low 
k, the dimension of this matrix is considerably smaller 
compared to that of an operator acting on the complete 
Hilbert space. Moreover, the size of the matrix depends 
only on the level k rather than the number of spins N. 

We next add a multiple of the total spin to and 
define the operators 

V^=C^ N)T^, (B2) 

where n is an integer. The operator is also closed in 
the subspace of states of level k if k < Ny since does 
not mix states of different levels. For certain values of 
n, we managed to find states constructed from current 
operators that are annihilated by the three operators 
for a G {x^y^z}. These states are then ground states of 
the Hamiltonian 


= (B3) 

where the index a is summed over. Note that the Hamil¬ 
tonian Hn is positive semi-definite and SU(2) invariant. 
Hn is nonlocal and contains terms with up to four-body 
interactions since has terms linear and quadratic in 
spin operators. 


since and commute. The operator + (I — 

N)T^ has the advantage that its matrix entries in terms 
of the states at level k do not depend N and n [cf. 
Eq. (C8) in Appendix [C]]. In our calculations, we found it 
technically easier to first determine the subspace of states 
on which and + (I — A)T^ commute and then 

look for states that are annihilated by for a suitable 
n within that subspace. 


TABLE 11. States constructed from current operators that 
are annihilated b y for a G {x,y, z} on a cylinder with 
Ny > k [cf. Eqs. ( |BI| and ( |B2| for the definition of P“]. For 
Ny sufficiently large {Ny > k), these states are ground states 
of the Hamiltonian Hn = 


k State 


'00 
Vi 

'03 + i£adei^\ 


2,1 


I a b 


Symmetric-traceless part of 03 ^ 
V’3,1.1 + |*eadeV’3.2 + |*£adeV’4,l + 


V>3 3 f + 


5 / a d d . 

3^421 \ 


'3 3 1 ^ T 4 2 1 3 T 4 2 1 

T 2 ^'^a(ie05^2 T T 




I 0 ? 


2 ‘^^aaer 6 ,l ' 4 

8 Symmetric-traceless part of 

iSade'lpl^l -h ^iSade'fpl^i^l “ §05,3 “ 06,2 “ 20?,^ 

9 Symmetric-traceless part of 0?,3,J 

I a d d 3 la dd _ I la d d 

2 V5,2,2 2V4,4,1 

,d d a 9jd ad _ 9ja d d _ 9jd a d 


9 ieadeipixij 


— dd _ 3 la 

2^4,3,2 2V'5 

+3^3,“ - iKlt - 

\ d d Q I a d d 


- ^ 7 /; 

2 V6,2,l 

d ( 


-ire,it - srit - fiSader,! - fi^adere, 

-^iSade'^il - ^iSadelpil “ 


Spin n 
0 any 

I I 

1 5 

2 3 

I 9 

1 13 

2 7 

3 5 

I 17 


We summarize our analytical results in Table \n\ The 
states with spin 2 and 3 appear as the symmetric- 
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traceless parts of states with 2 and 3 open indices, respec¬ 
tively. For a two-index state the symmetric-traceless 
part is defined as 

3{r^ + - 25ab4>'^^ (B7) 

and for a three-index state a^^ 

5 ((/>“'>'= -f -f (/>“'’ -f (/>'='’“ -f -f (/>“'='’) 

-f + </>"“)). (B8) 

Except for the levels 2 and 6, we find states and cor¬ 
responding parent Hamiltonians for all levels that were 
considered. Note that the singlet is a ground state of 
Hn for any value of n. For the additional ground states, 
we observe that the value of n tends to be larger at higher 
levels k. This means that the ground state space of the 
Hamiltonians with lower n contains states of a lower 
level in current operators. For example, we only find the 
ground states and 'ipi for Hi. Similarly, the only ap¬ 
pearing ground states of at levels k < 9 are 'ipo and 
the symmetric-traceless part of i. 

2. Ground-state degeneracies 

In the previous subsection, we explicitly constructed 
analytical ground states of the Hamiltonians with 
n G {1, 3, 5, 7, 9,13,17} in terms of linear combinations 
of states with levels k < 9. We now study the 

ground state spaces of the Hamiltonians numerically 
and provide evidence for n G {1,3,5} that the complete 
ground state space is spanned by the states given in Ta¬ 
ble HI 

TABLE III. Numerically determined ground state multiplets 
of the Hamiltonians Hn for n < 13 and an even number of 
spins N with N < 14. The second column indicates the mini¬ 
mal number of spins in the periodical direction for which 
the complete shown multiplet was observed in all system with 
< Ny < 14. For a lower number of spins in the y- 
direction, the observed ground state space is smaller. For 
even n, we only find a singlet ground state. 


n 


Ground state multiplet 

1 

2 

0© 1 

3 

4 

0® 2 

5 

6 

0© 1 ©3 

7 

8 

0©2©4 

9 

10 

0©1©3®5 

11 

12 

0®2©4®6 

13 

14 

0©1©3®5®7 


By an exact diagonalization, we numerically deter¬ 
mined the ground state multiplets of the Hamiltonians 
Hn for n < 13 and systems with N = N^Ny <14 and 
N even. Our results are summarized in Table imi We 
observe that states with spin s occur in the ground state 
spaces only in systems with Ny > 2s. Furthermore, we 
find that the ground state degeneracy does not grow any¬ 
more if Ny reaches a certain value N^^^. This statement 
is most conclusive for the lower values n, where N^^^ 
is smaller and we are thus able to probe more systems 
with Ny > N^^^. For n G {1,3,5}, this implies that all 
ground states are given by the corresponding states of 

Table m 

Finally, let us formulate a conjecture about the struc¬ 
ture of the states annihilated by which are ground 
states of Hn- Our analytical results are consistent with 
the following rule: For each spin sector s G {1,2,...}, 
there is a series of states at levels /c = 5^ + 2sj with 
j G {0,1,2,...}. These states are annihilated by 
with n = 25 — 1 + 4jf. As one can show by induction, the 
second rule implies that the ground state space of Hn 
with n = 25 — 1 contains the multiplet 

1 0 3 0 • • • 0 5, if 5 is odd, 

2 0 4 0 • • • 0 5, if 5 is even. 

The numerical results of Table HTTl are consistent with this 
multiplet structure and thus support the conjecture that 
the values of n are given by n = 25 — 1 0 4jf. 


Appendix C: Action of on states built from 
current operators 

Our starting point is the decoupling equation for the 
states derived in Ref.l33l This equation describes 

the action of the operator 

^i= E (Cl) 

on the states 'ipi^'/Pi^ and follows from the CFT null field 

{K^UJ^Mizi) with {K)i = 5ab-ieabct1. (C2) 

[The definition of Cf used here differs from that of Ref. [33] 
by a factor of 2/3.] 
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The operator 


N 


c- = J2ct 


(C3) 


i=l 

was used in Sec. to construct parent Hamiltonians for states built from current operators. 

The decoupling equation reads 

9=1 

+ 2{K^)i Y. (C4) 

Si,...,SAr g=2 n=0 


where 


^s(z) = . . .0SAr(^Ar). 


(C5) 


The decoupling equation for the states , where all mode numbers are one, is enough to describe the action of 

Cf on states with general mode numbers • Using the Kac-Moody algebra of Eq. the latter can be expressed 

in terms of the states ^ by repeated application of 




On the cylinder, we have 


N 


= 0, if n mod Ny ^ 0. 


i=l 


Summing over i in Eq. (C4), we therefore obtain for k < Ny 


(C6) 


(C7) 




q=l 


k q—1 


E EE ^a”..ai('Sl5 • • • ) SAr)|Sl, • • • , Sn), 

Si,...,S7v 9=2 n=0 


with 


Glz.a^ (si, ■■■,sn) = 2{M^){jx_,r\... rr • • • ^-\)(o)) 

- 25a,a{M^){Jln-lJ-\ • • • • • • ^-\)(0)). 

The first two terms on the right hand side of Eq. ( |C8[ ) are of order k in current operators since 

9=1 

In the remaining terms, the modes and can be commuted to the right since J^|0) = 0 for n > 0: 

= Y^saaAJ-V • • • j-TjLiJ-V • • • ^-\)(o)io), 


(C8) 


(C9) 


(CIO) 


(Cll) 


and similarly for U .. J^\)(0)|0). Iterating this step, the current operator modes with a positive mode number 

can be eliminated. The resulting terms only have negative mode numbers and are all of order k in current operators. 
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